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(N ; Abstract 

■ We consider the singlet scalar model of dark matter and study the expected 
Q\ , antiproton and positron signals from dark matter annihilations. The regions 
^p. \ of the viable parameter space of the model that are excluded by present data 

■ are determined, as well as those regions that will be probed by the forth- 
^ . coming experiment AMS-02. In all cases, different propagation models are 

H ! investigated, and the possible enhancement due to dark matter substructures 

is analyzed. We find that the antiproton signal is more easily detectable than 
the positron one over the whole parameter space. For a typical propaga- 
tion model and without any boost factor, AMS-02 will be able to probe -via 
antiprotons- the singlet model of dark matter up to masses of 600 GeV. An- 
tiprotons constitute, therefore, a promising signal to constraint or detect the 
singlet scalar model. 
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1 Introduction 



One of the simplest extensions of the Standard Model that can explain the dark 
matter is the addition of a real scalar singlet and an unbroken Z2 symmetry under 
which the singlet is odd while all other fields are even. Such a model, known as the 
singlet scalar model of dark matter, has been studied several times in the literature 
[D 121 El IH [6] . Predictivity is undoubtedly its most salient feature. In contrast with 
other common scenarios that explain the dark matter, the singlet model contains 
only one additional field, the singlet scalar, and two new parameters: the singlet 
mass and the coupling between the singlet and the higgs boson -the only standard 
model field that couples to it. The singlet relic density as well as its direct and 
indirect detections rates depend additionally on the higgs mass, which however 
cannot vary in a wide range. After imposing the dark matter constraint, the viable 
parameter space gets reduced simply to the singlet mass and the higgs mass. The 
implications of the model can then be studied as a function of the singlet mass for 
a few representative values of the higgs mass. The singlet scalar model, therefore, 
provides a compelling and predictive scenario to explain the dark matter. 

Recently, it was shown [6] that the singlet may have the right relic density to 
explain the observed dark matter abundance, that its direct detection cross section is 
large enough to be probed by present and planned experiments, and that the gamma 
rays from the annihilation of singlet scalar dark matter will likely be observed by the 
Fermi satellite [TJ. To complete the analysis of this model, in this paper we study the 
indirect detection of singlet scalar dark matter through positrons and antiprotons. 

The indirect detection of dark matter is one of the promising avenues towards 
its identification. Among them, antimatter searches in cosmic rays play an essential 
role. Present experiments such as Pamela [Sj and Fermi [7] are already measuring 
the cosmic ray spectrum with unprecedented precision. In fact they have already 
challenged our knowledge of the positron spectrum. Future experiments such as 
AMS-02 [9], which is scheduled to be launched next year, will measure the antiproton 
and the positron spectrum in a wider energy range and with significantly better 
statistics. Hence, the chances of observing an exotic component from dark matter 
annihilation in the positron or the antiproton spectrum are now higher than ever. 
It is therefore crucial to study the expected antiproton and positron fluxes in well- 
motivated scenarios for physics beyond the standard model that account for the 
dark matter of the Universe. 

The antiproton and positron fluxes from singlet annihilation not only depend 
on the dark matter model but also on a number of astrophysical factors that affect 
the production and propagation of such particles throughout the Galaxy. To take 
these effects into account, we consider different propagation models for positrons and 
antiprotons, and we study the possible role of substructures in the dark matter halo. 
Our goals are to find out whether any region of the parameter space is already ruled 
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out by present data from Pamela, and to identify those regions that are within the 
sensitivity of the AMS-02 experiment. In addition, by comparing the corresponding 
regions for positrons and antiprotons, we expect to determine which of these two 
signals offers better perspectives to probe the singlet scalar model. 

In the next section the main features of the singlet scalar model of dark mat- 
ter will be briefly reviewed. Then, in section [31 we present a detailed analysis of 
the antiproton signal in this model. In particular, we determine the excluded and 
detectable regions in the plane (mg, (crv)) for different sets of parameters account- 
ing for propagation and substructure effects. In section 0] we present an analogous 
analysis for the positron channel. A brief discussion of our main results is given in 
section \5[ which is followed by our conclusions. 

2 The singlet scalar model of dark matter 

The Lagrangian that describes the singlet scalar model of dark matter is 



where Csm denotes the Standard Model Lagrangian, H is the higgs doublet, and S 
is the singlet scalar field. This is the most general renormalizable Lagrangian that is 
compatible with the SU(3) x SU(2) x U(l) gauge invariance and with the symmetry 
S — > —S. Notice, from the Lagrangian, that the higgs boson is the only standard 
model field that directly couples to the singlet. This extension of the standard 
model contains, therefore, two new phenomenologically relevant parameters: mo 
and A. Instead of mo, it is useful to consider the physical mass of the singlet 



as a free parameter of the model. The other free parameter is A, which determines the 
trilinear (S 2 h) and quartic (S 2 h 2 ) coupling between the singlet and the physical higgs 
boson. In addition to mg and A, the dependence on the higgs mass -a SM parameter- 
should also be also taken into account, as it affects the dark matter phenomenology. 
To accurately compute the singlet relic density we use the micrOMEGAs package 
|10j . which can calculate the relic density in a generic dark matter model. 

Singlets can annihilate through s-channel higgs boson exchange into a variety 
of final states: //, W + W~, Z°Z°, and hh. Additionally, they can also annihilate 
into hh either directly or through singlet exchange. Depending on the singlet mass, 
two annihilation regions can be clearly distinguished. A light singlet, mj < mw, 
will annihilate mainly into the bb final state, with subdominant contributions from 
other light fermions. To obtain the correct relic density, such light singlets require 
a relatively large value of A (> 0.1), and are consequently constrained by direct 



C = C SM + ^d,Sd»S - ^S 2 - ^S 4 - \S 2 H^H , 
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Figure 1: Annihilation branching fractions as a function of the dark matter mass 
for the singlet model. The higgs mass was set to 120 GeV while the value of A was 
obtained by imposing the dark matter constraint. 



detection experiments. Singlets masses below 50 GeV, for instance, are already 
ruled out by current measurements [6j. A heavy singlet, ms > rriw, on the other 
hand, will annihilate mainly into W + W~, with additional contributions from Z°Z°, 
hh and tt. In this region the required value of A is typically smaller (~ 0.01) and 
present direct detection constrains are ineffective. Figure [T] displays the annihilation 
branching fractions as a function of the singlet mass for = 120 GeV. The sharp 
contrast between the light singlet region and the heavy one is clearly observed. 

Being a scalar field, the singlet annihilation rate today is not suppressed with 
respect to that of the early Universe (s-wave annihilation). Moreover, throughout 
most of the viable parameter space (av) is constant and equal to the typical annihi- 
lation rate, (av) ~ 2.3 x 10 _26 cm 3 s _1 [BJ. The only two places where (av) is smaller 
are at the W threshold (mg = mw — e) and at the higgs resonance (m^ ~ 2ms). 
In fact, the main effect of the higgs mass is to determine the position where this 
resonance lies. In the following analysis we will set rrih to 120 GeV and consider 
only the remaining dependence on the singlet mass. For the sake of definiteness, 
singlet masses below 600 GeV will be considered. 
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3 Antiprotons 



Antiprotons can be produced either by dark matter annihilations or through different 
kinds of astrophysical mechanisms. After being produced they propagate in the 
Galaxy and reach the Earth, where they can be detected as exotic components in 
cosmic rays. The PAMELA experiment reported recently [H] the measurement of 
the antiproton to proton flux ratio up to 100 GeV, which is in agreement with that 
predicted by the conventional background model. In the near future, the AMS-02 
experiment, to be launched next year, will measure the antiproton flux with even 
higher precision, increasing the odds of finding an additional component due to dark 
matter annihilations. 



3.1 Propagation 

Once produced, antiprotons propagate throughout the galaxy undergoing several 
processes: 

• They scatter on irregularities of the galactic magnetic field -Alfven waves. 
These scatterings constitute, in fact, a space diffusion process with a diffusion 
coefficient given by 

K(E kin ) = K (3 p (-^) S , (3) 

where E kin is the antiproton kinetic energy, p p = (E kin + 2m p E kin ) 1 ^ 2 is the 

antiproton momentum, and 8 P — 1 — — -f — w ■ K and 5 are free pa- 
rameters of the propagation model that are constrained by a combination of 
theoretical predictions and astrophysical data. 

• They undergo a second order Fermi mechanism reacceleration, due to the 
motion of the scattering centers with a velocity of V a ~ (20 — 100) km/sec. 
This reacceleration process is described by the coefficient 

K EE = ~V 2 ^. (4) 
EE 9 a K(E) 1 1 

• They lose energy either adiabatically, or through Coulomb scattering, or by 
ionizing the interstellar (IS) medium. The total energy loss rate is denoted by 
b and depends on the antiproton energy. 

• They are wiped away from the galactic disk through convection, with a velocity 
V c ~ (5 — 15) km/sec. In the following, this velocity will be taken to be 
completely vertical to the disk: V c = V c sign(z). 
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• They can annihilate upon scattering on the IS medium. In this study, we shall 
consider the two primary components of the medium, namely Hydrogen and 
Helium. The annihilation cross-sections for p — H and p — He scattering have 
been taken from [12], where the well-known Tan&Ng parametrization [13] is 
used. 

The propagation of antiprotons in the interstellar medium, taking into account all 
of the above processes, is described by a diffusion-convection equation of the form 

d z (V c ^) - KV^j + d E [b{E)i) - K EE (E)d E ^} =q, (5) 

where we denote by ip = dn/dE the energy density of the antiprotons, and by q 
the source term -see equation (|U]). Actually, this equation is rather generic and can 
be applied to other particles species propagating in the galaxy. Only the values 
of the different parameters, but not the equation itself, will vary depending on the 
propagating particle. We will see in section @J for instance, that this same equation 
describes the propagation of positrons in the interstellar medium. Let us now see 
how to solve this equation for antiprotons. 

3.2 The primary flux 

To solve equation (jSJ), we shall use the method proposed, for example, in [Ti t [TBI, [T5] . 
The idea of this method is to adopt a simpler version of equation (j5j) for which 
the Green function can be calculated analytically. To achieve such simplification, 
certain processes contributing to the final antiproton spectrum have to be neglected. 
Specifically, all energy redistributions in the initial (injection) spectrum -energy 
losses, reacceleration, as well as "tertiary" contributions (i.e. contributions from 
secondary antiprotons produced upon inelastic scattering with the IS medium)- are 
ignored. Whether these redistributions are important or not depends mainly on the 
antiproton energy. For GeV energies, the results may deviate up to 50% from those 
obtained with the (more complete) Bessel function treatment -in [19] a comparison 
between the two methods can be found (see figure 2). But for energies around 
10 GeV, the accuracy of the method improves dramatically, yielding essentially 
indistinguishable results at slightly higher energies. Since the p energy region we 
shall consider begins at 10 GeV, we can safely use this simplified approach. Apart 
from its clearer physical interpretation, the main advantage of this method is that 
astrophysical effects can be separated from particle physics ones. As a result, it is 
a well-suited method for scans in parameter spaces, as those we are going to carry 
out in the following. 

If we denote by r^ nn = ^2 1SM n lSM vap™ SM the destruction rate of antiprotons in 
the ISM, where ISM = H, He, and implementing the aforementioned simplifications, 
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the transport equation for a point source (which actually defines the propagator G) 

G = 5(f-r y ). (6) 
The antiproton propagator can then be written as 



is: 

d 



-KV + V c — + 2hT tot 5(z) 



G *^ Z) = ^^Y, c n lK °( r VW^JMk n L)sm(k n (L-z)), (7) 



2nKL 

n=0 

where K is a modified Bessel function of the second kind and 

^ sin(fc n L) cos(fc n L) 
k n L 

k v = V C /(2K), (8) 
k d = 2hV^/K + 2k v , 

with h = 100 pc being the half thickness of the galactic disc, and L being the 
half-thickness of the diffusive zone. k n is obtained as the solution of the equation 

mi — k n L — arctan(2fc n /fcrf) =0, iiGN. 

Then, in order to compute the flux expected on earth, we should convolute the Green 
function (J7J) with the source distribution g(r*, E). For dark matter annihilations in 
the galactic halo, the source term is given by 

^>4(f) 2 x:Kf)> 

where the index i runs over all possible annihilation final states. The decay of SM 
particles into antiprotons can be calculated with programs like PYTHIA [22J. In 
the singlet scalar model of dark matter, antiprotons may orig inate in W ± , Z°, h, 
and t decays. Regarding the distribution of dark matter in the Galaxy, p(x), we 
assume a NFW profile with a local density of 0.3 GeVcm 3 . The final expression for 
the antiproton flux on the Earth takes the form 

- S¥ (€r) 2 >"°> L (Mi G ^> (10) 

where none of the integrated quantities depends on the antiproton energy. This 
feature demonstrates one of the virtues of the Green function method applied to an- 
tiprotons: The integral in equation (1101) . which we compute using a VEGAS Monte- 
Carlo algorithm, needs to be calculated only once for each value of the injection 
energy, for it is the same as the detection energy. 
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L (kpc) 


K (kpc 2 /Myr) 


S 


V c (km/s) 


MIN 


1 


0.0016 


0.85 


13.5 


MED 


4 


0.0112 


0.70 


12.0 


MAX 


15 


0.0765 


0.46 


5.0 



Table 1: Values of propagation parameters widely used in the literature and providing minimal 
and maximal antiproton fluxes, or consitute the best fit to the B/C data. 

Regarding the propagation parameters L,K ,6, and V c , we take their values 
from the well-established MIN, MAX and MED models -see table [TJ The former 
two models correspond to the minimal and maximal antiprotons fluxes that are 
compatible with the B/C data. The MED model, on the other hand, correspond to 
the parameters that best fit the B/C data. 



3.3 Influence of substructures 

N-body simulations reveal that galactic halos are not completely smooth; they also 
contain a significant number of substructures (clumps). Such substructures have 
been studied repeatedly in the literature as a possible way to enhance the dark 
matter annihilation rate. It has been claimed though, that in realistic scenarios it 
is rather improbable to expect large enhancements from dark matter clumps [T4] . 
Here, we will study rather qualitatively the possible effect of substructures on the 
antiproton (and positron) signal. 

To that end, we closely follow the approach outlined in [15J. The enhancement 
due to dark matter substructure is described by an energy-dependent function known 
as the boost factor B. Because the distribution of dark matter clumps in the Galaxy 
is unknown, B cannot be computed from first principles; it can only be studied 
from a statistical point of view. For simplicity, we will limit our discussion to the 
effective boost factor, B e ff, defined as the average value of B over a large number 
of realizations of the Galactic halo. It must be kept in mind, however, that in some 
exceptional cases -for example when there is a large dark matter clump very close 
to the Earth- -B e // can deviate significantly from B. Since the probability of such 
an event is quite small [17], we will not consider those cases in this work. 

The energy-dependent effective boost factor, -B e //, can be written, under certain 
assumptions (see [IS]), as 

B CS EE M = (1 _ ff + f Bc ^ , (11) 

where ((f)) is the average flux coming from the clumpy DM distribution, sm is the 
flux that we would expect if the whole halo were smooth, / is the fraction of DM in 
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clumps, and B c is the boost factor (assumed constant with energy) for an individual 
clump. In this study, this constant boost factor is supposed to be universal for all 
clumps. Finally, I n =i,2 are given by 



The effective boost factor, then, depends on / and B c . When invoking dumpiness, 
we will follow [15] and use / = 0.2 as a representative value (see e.g. [201 [21]). Re- 
garding the constant boost factor, B c , it could vary from just a few up to two orders 
of magnitude [HI [TBI El]. We will use B c = 3, 10, 100, which give rise to effective 
boost factors in the approximate ranges (1,2), (3,5) and (10,40) respectively. This 
last range roughly coincides with the upper limit for the boost factor found in [14] , 
for the case of a NFW smooth halo and clumps with a Moore et al internal profile. 

3.4 Astrophysical backgrounds 

So far, our main concern has been the antiproton signal from dark matter annihila- 
tions. But, to determine whether a model is detectable or whether it is compatible 
with present data, we also need to know the antiproton background. The most 
well-known treatment of the astrophysical antiproton background is that of Strong 
and Moskalenko [T2l . which has recently received additional experimental support. 
In fact, the antiproton data [11] from the PAMELA experiment is perfectly compat- 
ible with their predictions. We thus have a relatively simple antiproton background 
model that is able to explain the data. 

In the following, when trying to examine whether the model is excluded by 
the PAMELA data, we shall be employing the theoretical background predictions, 
normalizing them to the low-energy PAMELA data. More specifically, we take 
advantage of the fact that predictions demonstrate a relatively constant background 
with respect to the propagation model, and use the simple parametrization (properly 
normalized) presented in [19J. When studying the detectability of the model in 
AMS-02, on the other hand, we shall use the PAMELA measurements themselves 
as background. More concretely, we borrow the fit performed by Cirelli et al in [23J , 
which provides a sufficiently good fit to theoretical predictions but also to the recent 
PAMELA data. We pay special attention at reproducing the good background 
normalization at low energies, so as to stay as close as possible to the PAMELA 
measurements. 

Finally, let us point out that since the antiproton background does not vary 
significantly with the propagation parameters [T9] , it is reasonable to consider a 
global background, irrespective of the propagation model. We will do so in the 
following. 
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Figure 2: Regions of the parameter space that are excluded by the recent antiproton data from 
the PAMELA experiment. The area above the MIN, MED, and MAX lines is excluded for the 
given propagation model. The solid (black) line shows the viable parameter space of the singlet 
scalar model of dark matter. 



3.5 Results 

In this section we present our results for antiprotons. First, we combine the predicted 
flux and the expected background to determine the regions in the plane (ms, (av)) 
that are already excluded by present antiproton data from PAMELA. Then, we 
determine the regions that are detectable by AMS-02. 

In the singlet scalar model of dark matter, the singlet mass and the higgs mass 
are the only free parameters left after imposing the relic density constraint. Since 
we take rrih = 120 GeV throughout this work, the annihilation branching ratios and 
the total annihilation rate, (av) sing i et , depend only on ms- To obtain the excluded 
regions in the plane (ms, (&v)), we first compute, for a model with the same branch- 
ing ratios as the singlet scalar model (see figured]), the value of (av) required to 
exclude the model, (&v) excl . By comparing (&v) singlet with (cv) exd we can then de- 
termine whether the model is excluded or not at a given singlet mass. An analogous 
procedure is followed to determine the detectable regions. 

The condition we impose for exclusion is that the sum of the signal prediction 
and the (properly normalized) background prediction exceed the measurement by 
the PAMELA collaboration by more than 3a. Figure [2] shows the excluded region 
for the MIN, MED, and MAX propagation models. The area above the lines is 
already excluded by present data. Also shown as a solid (black) line is the region 
along which the relic density constraint is satisfied -the viable parameter space of 
the model. Because that line is entirely below the exclusion lines, we conclude that, 
in the absence of substructure enhancements, current antiproton data do not yet 
constrain the singlet scalar model of dark matter. This situation may change soon, 
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Figure 3: Regions of the parameter space that are within the sensitivity of the AMS-2 experiment. 
The area above the MIN, MED, and MAX lines is detectable by AMS-02. The solid (black) line 
shows the prediction of the singlet model. Notice that for MED and MAX, essentially the whole 
parameter space is detectable. 

however, for the MAX exclusion line lies just above the prediction of the model for 
masses around 100 GeV. Future data from PAMELA may start constraining that 
region of the parameter space. 

To assess the sensitivity of AMS-02 to an antiproton signal from singlet dark 
matter annihilations, we compute the annihilation rate needed to distinguish the 
signal from the background at the 95% confidence level. To do so, we first calculate 
the x 2 as 



where <p tot is the total antiproton flux, <j) bk9 is the background flux, N is the number 
of energy bins considered, A is the geometrical acceptance of the experiment, and 
T is the data acquisition time. AMS-02 is expected to take data for three years 
and features an antiproton geometrical acceptance of 330 cm 2 sr [21]. We consider 
20 energy bins evenly distributed in logarithmic scale between 10 and 300 GeV. A 
95% confidence level corresponds then to x 2 > 31. 

The annihilation rates needed for a detection at AMS-02 are shown in figure 
[3] for the three propagation models. The singlet model is detectable whenever its 
annihilation cross-section (solid line) is larger than the exclusion lines at AMS-02 
for a given propagation model. We see that for the MIN propagation model only a 
small region around 60 GeV is within the AMS-02 reach. For the MED and MAX 
propagation models, on the other hand, essentially the whole parameter space is 
detectable at AMS-02. Only in the higgs resonance region, 2ms ~ ^h, and at the 
iy-threshold does the model's annihilation rate fall below the exclusion lines. Hence, 
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for the MED and MAX models, a singlet scalar with a mass below 600 GeV should 
be easily detected at AMS-02. 

We now consider the possible effect of dark matter substructures on the exclusion 
and the detectable regions. Figure H] shows the exclusion lines for the MIN (upper 
figure), MED (middle figure), and MAX (lower figure) propagation models for the 
different values of B c we investigate. As before, the solid (black) line shows the 
prediction of the singlet model. Notice that the MIN propagation model is nowhere 
excluded, not even for the most optimistic enhancement. For the MED propagation 
model, B c = 100 is excluded for singlet masses below 360 GeV. The other values of 
B c are all compatible with present data. For the MAX propagation model, B c = 100 
is ruled out in the whole mass range while B c = 10 is at the limits of exclusion for 
masses around 100 GeV. Hence, we see that, once the possible enhancement due 
to dark matter substructures is included, current antiproton data already constrain 
some combinations of astrophysical and particle physics parameters. 

The prospects for the detection of an antiproton signal at AMS-02 were already 
good without any additional boost from dark matter clumps, as we saw in figure 
El The additional enhancement from dark matter substructures only improves such 
prospects, as clearly seen in figure For the MIN model, for instance, the detectable 
region increases from a small region around mj ~ 60 GeV without boost factor to 
m s ^ 200 GeV for a moderate boost factor, covering the whole mass range for 
B c = 100. 

To summarize, we have seen that the antiproton signal from singlet annihilation 
is a promising way to indirectly detect singlet scalar dark matter. Current antiproton 
data from PAMELA already constraint some scenarios and future data from AMS-02 
is likely to reveal a significant excess over the expected background. 

4 Positrons 

We now discuss the positron signal from singlet dark matter annihilation, carrying 
out an analysis similar to that presented for antiprotons in the previous section. 

4.1 Propagation 

There have been many treatments on the propagation of positrons throughout the 
galactic medium. In this paper, we adopt the two-zone diffusion model and solution 
described in [15]. For completeness, we review here the main points of this approach. 

Being a diffusive process, positron propagation is governed by the same gen- 
eral equation, that describes antiproton propagation. The approximations and 
the parameters that were used for antiprotons, however, are not the same as for 
positrons. Indeed, contrary to the p case, where energy redistribution processes be- 
come inefficient above a few GeV, energy loss -through either synchrotron radiation 
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Figure 4: Regions excluded by the antiproton data from PAMELA including the possible effect 
of substructures in the DM halo. From top to bottom the figures correspond to the MIN, MED, 
and MAX propagation models. The solid (black) line shows the prediction of the singlet model. 
The area above the lines is excluded for the corresponding parameter values. 
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Figure 5: Detectable regions at AMS-02 including the possible effect of substructures. From top 
to bottom the figures correspond to the MIN, MED, and MAX propagation models. The solid 
(black) line shows the prediction of the singlet model. The area above the lines is detectable for 
the corresponding parameter values. 
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L (kpc) 


K (kpc 2 /Myr) 


a 


MIN 


1 


0.00595 


0.55 


MED 


4 


0.0112 


0.70 


MAX 


15 


0.0765 


0.46 



Table 2: Values of propagation parameters widely used in the literature and roughly providing 
minimal and maximal positron fluxes, or constitute the best fit to the B/C data. 



or inverse Compton scattering on stellar light and CMB photons- is the main pro- 
cess involved in positron propagation. These processes lead to an energy loss rate 
that can be written as 

b(E) =€i' (i4) 

where E is the positron energy, Eq is a reference energy (which we take to be 1 
GeV) and te = 10 16 s is the characteristic energy-loss time. 

The other important difference with respect to antiprotons is that for positron 
propagation the effect of the galactic convective wind, as well as reacceleration and 
annihilation processes can all be neglected. After these simplifications are taken into 
account, we are left with the following equation 

d t iP - V [K(x, E)VtP) - d E [b(E)iP] = q(x, E) , (15) 

where K is the space diffusion coefficient -steady state is assumed. This coefficient 
is taken to be constant in space but depends on the energy as 

K(E) = K (JL) . (16) 

Here the diffusion constant, K , and the spectral index, a, are propagation parame- 
ters. Then, the master equation for positron propagation gets simplified to its final 
form 



K e a A^ + -^—^j + g = 0, (17) 

where e = E/E . This is the equation we solve to calculate the effects of positron 
propagation on a signal produced at some point in the galaxy. 



As in the case of antiprotons, a crucial factor in the treatment of positron prop- 
agation is the adopted propagation model. In this case there are mainly 3 relevant 
parameters, namely L, K and a; that is, the half-thickness of the cylindrical dif- 
fusive zone, the diffusion constant and the spectral index, respectively. For their 
values we use the three models described in table [2j 
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4.2 The primary flux 

The resulting positron flux from DM annihilations can be written as (see [25J for 
details) 

where the detection and the production energy are denoted respectively by E and 
E s . f(E s ) is the production spectrum for positrons, f(E s ) = '^2 i dNl + /dE s , with i 
running over all possible annihilation channels. The diffusion length, A^, is defined 
by 

(pa- 1 fpa—1 \ 
■ < 19 > 

The so-called halo function, /, contains all the dependence on astrophysical factors. 
It is given as 

I(X D ) = [ d 3 x s G(^ E -> x s , E s ) ( ^A) • ( 20 ) 
Jdz \p{xq)J 

The modified Green function G is in its turn defined by 



G = — - — e -(*Q-zs) 2 /*Ko?y (21) 



with V depending on the value of the characteristic parameter £ = When 
C > 1 -when the diffusion time is small- boundary conditions can be ignored and 
the propagation equation can be treated as a 1-D Schroedinger equation. In that 
case 

When £ is small this approximation no longer holds but we can express V as 



V = 

n=l 



oo 

J2 J [ e ~ Xnf( f>n(Zs)<f>n(z) + e _At Vn(2.KW] (23) 



sin(A; n (L - \z\)) , fc n = ( n - ^ ) ^ (24) 



where 



2/ L 

= sin(A£(L - *)) , *4=™f- (25) 

We now have all the necessary ingredients for the computation of the expected 
positron flux from singlet scalar dark matter annihilation. Next we discuss the 
positron background. 
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4.3 Astrophysical Backgrounds 



In the conventional background model, positrons are produced in the interactions 
between cosmic-ray nuclei and the interstellar medium. The expected positron flux 
on earth can be written as 126 



4 5 E - 7 

§ C X V = ^ GeV^cnT V^r" 1 (26) 

1 + 650£ 2 - 3 + 1500E 4 - 2 v ; 

where E is given in GeV. In contrast to the antiproton case, this conventional 
background is not compatible with the recent data from PAMELA [27J and Fermi- 
LAT [28] . Even after taking into account the possible uncertainties due to cosmic ray 
propagation, the data reveals a clear excess over this background at high energies, 
E > 10 GeV. Hence, a new source of high energy positrons is necessary to explain 
the data. 

Two main interpretations of this positron excess have been considered in the 
literature. One is the annihilation of dark matter particles. Though viable, such 
models are rather contrived. They require multi-TeV masses, annihilation into lep- 
tonic final states, and an annihilation rate several orders of magnitude larger than 
that expected for a thermal relic [29J. Moreover, they are tightly constrained by 
experimental data -gamma rays, radio emission, CMB, BBN [30]. Given that, in 
addition, the singlet model cannot explain the positron excess (at least in the en- 
ergy region we examine and for reasonable astrophysical assumptions), we will not 
consider this dark matter interpretation any further. 

The other main interpretation for the observed positron excess is that it is gen- 
erated by astrophysical sources, such as pulsars [21] • If that is the case, such astro- 
physical positrons constitute an additional background to the positrons from dark 
matter annihilation. And they must be taken into account when studying the sen- 
sitivity of future experiments to a dark matter positron signal, as pointed out in 
[32J. Here, we will closely follow their approach. That is, we use as background the 
positron flux obtained from the pulsar interpretation of PAMELA and Fermi-LAT 
data [33] -see figure 2 in [32 j . That background is compatible with the absolute 
positron flux derived in [3] 



4.4 Results 

For positrons we follow a procedure analogous to that for antiprotons. We first find 
the regions that are already excluded by present data and then determine those 
that will be detectable in AMS-02. When using equation f[T3|) . it must be taken into 
account that the geometrical acceptance of AMS-02 for positrons is 420cm 2 sr [24J. 

Figure [6] displays the regions that are already excluded for the MIN, MED, and 
MAX propagation models. The solid (black) line shows, instead, the prediction of 
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Figure 6: Regions of the parameter space that are excluded by the recent positron data from 
the PAMELA experiment. The area above the lines is excluded for the corresponding propagation 
model. Notice that no region of the viable parameter space is currently ruled out. 
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Figure 7: Regions of the parameter space that give a positron signal within the sensitivity of the 
AMS-02 experiment. The lines corresponding to the MIN, MED and MAX propagation models 
are shown. They must be compared to the actual prediction (solid line) of the singlet scalar model. 



the singlet scalar model. Notice that the variation due to the propagation model is 
very small, certainly much smaller than for antiprotons. From the figure we also see 
that the exclusion lines all lie well above the model's prediction. We checked that a 
similar result is found even in the optimistic scenario for substructure enhancement. 
At present, henceforth, no region of the viable parameter space is even close to being 
excluded via the positron signal. 

The AMS-02 detectable regions are shown in figure [7] for the three propagation 
models. Again, not much difference is observed between the MIN, MED and MAX 
lines. Moreover, only a small region around mg ~ 200 GeV is within the AMS-02 
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reach. Thus, without additional boost factors, singlet masses above 200 GeV do not 
produce a detectable positron signal at AMS-02. 

Once the effect of dark matter substructures is included, the detectable region 
extends to higher singlet masses, as illustrated in figure For the optimistic sce- 
nario, B c = 100, the whole mass range becomes detectable independently of the 
propagation model. For the moderate scenario, B c = 10, the AMS-02 reach extends 
to singlet masses of 100 GeV and 200 GeV for the MIN and MAX models respec- 
tively. Thus, a substructure enhancement is required to probe the singlet model 
beyond m s ~ 200 GeV at AMS-02. 

5 Discussion 

Our results clearly indicate that in the singlet scalar model the indirect detection of 
dark matter via antimatter is feasible, particularly through the antiproton channel. 
The positron channel also gives rise to an observable signal at AMS-02, but only for 
light singlets. The antiproton signal, on the other hand, should be visible over most 
of the viable parameter space. 

To identify the singlet scalar model as the correct model of dark matter, one 
would need to measure rrih and ms, the only two parameters of the model, and to 
confirm some of the model predictions. Such a task would certainly require input 
from different experiments: accelerator searches, direct detection experiments, and 
indirect searches through several channels. Searches at the LHC, for instance, could 
provide a measurement of the higgs mass, but they are not expected to say much 
about the singlet mass. It is mainly through direct and indirect detection searches 
that the singlet mass could be determined. 

In [BJ, it was shown that the singlet direct detection rate is within the reach of 
current and future experiments and that the gamma ray signal from singlet dark 
matter annihilation is likely to be observed by the Fermi satellite. In this paper, we 
have found that also the antiproton signal is observable at AMS-02. By combining 
the signals from these three experiments the singlet mass, ms, might be determined. 
Once m/i and ms are known, the singlet model of dark matter is completely specified, 
at least within the standard cosmological scenario. 

6 Conclusions 

The singlet scalar model provides a simple and predictive scenario to account for the 
dark matter. In this paper, we studied in detail the antiproton and positron signals 
expected from the annihilation of singlet scalar dark matter. Particular attention 
was paid to the role of the propagation model and to a possible enhancement in the 
annihilation rate due to dark matter substructures in the Galaxy. We used the recent 
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Figure 8: Detectable regions for the MIN, MED, and MAX propagation models in the presence 
of dark matter substructures. The region above the lines is detectable for the corresponding 
parameters. The solid (black) line shows the prediction of the singlet scalar model. 
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measurements, from PAMELA and Fermi-LAT, of the positron and antiproton flux 
to find the regions in the plane (m s , (cry)) that are already excluded by the data. 
For antiprotons, an optimistic enhancement, B c = 100, due to dark matter clumps 
is essentially excluded for both the MED and the MAX propagation models in the 
entire mass range we considered. For positrons, on the contrary, no region is excluded 
independently of the propagation model or the possible effect of substructures. We 
also analyzed the perspectives for detecting antimatter signals in future experiments. 
In AMS-02, the positron signal, in the absence of additional boosts, was found to be 
detectable but only for light singlets. Antiprotons, on the other hand, offer excellent 
perspectives to be detected at AMS-02. Even without boosts, the whole mass range 
is within the sensitivity of AMS-02 for the MED and the MAX propagation models. 
For the MIN model, a moderate enhancement from substructures, B c = 10, would 
allow to probe the mass range below 200GeV. Antiprotons, therefore, seem to be a 
promising way to constrain or detect the singlet scalar model of dark matter. 
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